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Abstract 

The analogue of the Kolmogorov's four-fifths law is derived for two- 
dimensional, homogeneous, isotropic EMHD turbulence in the energy cascade 
inertial range. Direct numerical simulations for the freely decaying case show 
that this relation holds true for different values of the adimensional electron 
inertial length scale, de- The energy spectrum is found to be close to the 

expected Kolmogorov spectrum. 
PACS numbers: 47.27, 52.35 
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The statistical theory of three-dimensional fully developed hydrodynamic turbulence re- 
lies on one outstanding issue: the nonlinear transfer of energy from large to small scales 
It is therefore interesting to look for two-dimensional turbulent fluid dynamical systems shar- 
ing this same feature. Actually many of them exhibit a reversed energy flux, from the small 
scales to the larger ones, as is the case of 2D Navier-Stokes turbulence Hasegawa-Mima 
turbulence |§| or its geophysical counterpart, equivalent barotropic turbulence 

In this framework 2D electron-magnetohydrodynamic (EMHD) turbulence deserves spe- 
cial attention, beyond its modeling applications, since it has been shown to display, for the 
freely decaying case, a forward energy cascade a la Richardson-Kolmogorov ||10|| . 

In this Letter it is introduced a relation which is the counterpart of the Kolmogorov four- 
fifths law for homogeneous and isotropic 2D EMHD turbulence. Its content is compared with 
the results obtained by direct numerical simulations. 

EMHD equations are a fluid dynamical model for a cold electron plasma, moving in 
a uniform charge-neutralizing background of stationary ions. In recent years this model 
has received considerable interest for its relation to inertially confined plasma and to laser- 
plasma interactions, but the comparison with experimental results is limited by the fact that 
plasma which evolve according to EMHD equations are usually short-lived. 



The equations for the electron plasma in adimensional form are |Tl 



dl— + dlv ■ Vv = -Vp -E-VXB + Ug{-V)^''v . (1) 

The velocity field v, the magnetic field B, and the electric field E are allowed to have a 
nonzero component in any direction of three-dimensional space, meanwhile their functional 
dependence is restricted to the plane coordinates {-^ = 0). The equations have been adi- 
mensionalized with respect to the characteristic macroscopic length L, the typical magnetic 
field -Bq, the characteristic time r = mf.c/ {eB^d^^) and rfg = [meC^/(47re^nL^)]^/^ is the ratio 
of the inertial electron length scale to the integral scale L. The density of the number of 
electrons n is assumed to be uniform according to the incompressibility of the velocity field 
V-f = 0. The approximation made by considering motionless ions requests that the adimen- 



sional ion inertial length, di, must be larger than unity, thus limiting the range of admissible 
values for de = di{rae/m.i)'^^'^ to the interval d^. > 0.02. As long as we are dealing with a 
fluid description of plasma, all the lengthscales under consideration must largely exceed the 
Larmor radius. For stationary ions and negligible displacement current, the Ampere law 
becomes 

V = -V X S . (2) 

The last term in equation (P is a dissipative term which mimics the effects of electron 
viscosity (g = 1) or resistivity (g = 0). The total energy of the electron fluid 

E = ]^j {dlv^ + B^)d''x (3) 

is conserved by these equations in the ideal, non-coUisional case. 

In the spirit of Kolmogorov analysis one takes under consideration the spectral 

energy budget 

^^S{K) + Ii{K) = -V{K) (4) 

where £{K) = E{k)dk is the mean cumulative energy per unit mass contained at 
wavenumbers smaller than K, and E{k) is the energy spectrum. T{{K) is the energy flux 
(per unit mass) from wavenumbers k < K to larger wavenumbers, and "D^K) is the cumu- 
lative energy dissipation up to wavenumber K. Since the dissipation is localized to high 
wavenumbers k > K^, there is a range of wavenumbers Kq <^ K <^ where T>{K) ~ 0, 
and the energy flux Il{K) is determined by the inertial transfer of energy from the energy 
containing eddies at wavenumbers around Kq. When the large scale energy input due to 
the straining of energy containing eddies at scale Kq is equilibrated by dissipation taking 
place at small scales, one expects the energy flux through wavenumber K, Il{K), to be 
independent of i^' 0. Actually it must be remarked that, due to the energy decay, the 
flux approaches a constant value, equal to the total dissipation, only for very large K, and 
the crossover to the asymptotic behavior is very slow [|12|. In an analogous fashion as in 
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three-dimensional hydrodynamic turbulence, in the limit of vanishing viscosity, z/ ^ 0, the 
energy flux is expected to achieve a flnite positive limit, depending on the value of cie, ^ > 0, 

U{K) Ko<^K<^Kd. (5) 

This is a strong request which in 3D-hydrodynamics has experimental evidence; in the case 
under consideration it will be shown that numerical simulations provide reasonable support 
to this hypothesis. 

Assuming statistical homogeneity the energy flux, Il{K) can be expressed by means of 
physical space statistics by performing the Fourier transform of (Q). The result is 

U{K) = ^ I K^^^e{l) . (6) 



27r 

where Ji is the flrst-order Bessel function of the flrst kind, and the energy flux in the physical 
space is given by 

'^^^ = "^^^'^'''^^^ ■ '"^'^ + + ^^^^ ■ + '^^ 
where the subscript NL stands for the nonlinear contribution to the time derivative of the 
fields, as it can be extracted by the equations of motion and the brackets (. . .) express 
ensemble averages. Using ([l|) and (|), making repeatedly use of statistical homogeneity, of 
incompressibility of the velocity field and solenoidality of the magnetic field, one obtains the 
following relation 

9£A 4 (8) 
-({5v ■ 5B)5B,) - -{{5B ■ 5B)5vi)) 

with z = 1, 2 denotes the planar components. The expression (H) for the physical space 
energy fiux is the analogue of the Karman-Howarth-Monin relation |0] , and it involves only 
differences of dynamical fields as 5v = v{x -\- £) — v{x). To proceed further one assumes 
statistical isotropy and it is then possible to show that the physical space energy flux is 



in which appear the following third order structure functions 

Us{i) = {6v\\6BJBj_) , V^{i) = {6vs6v;,6v\\) , (10) 

Ws{e) = {6vs6B^s6Bii) , X3(£) = {6B^6B^6vii) , 

where the standard notation for longitudinal, 6v\\ = 6viii/i, and transverse differences, 
Sv± = eij6viij/£, has been used. 

As a consequence of hypothesis (|^), it can be shown by a saddle-point argument that 
the physical space flux, e{i), must behave as 

e{i) ~ e A < £ < 4 , (H) 

in the limit of vanishing viscosity, where the inertial range of lengthscales is now delimited 
by the "Taylor scale", A = {ugE/Dy/^i, where D is the energy dissipation, and the energy 
containing scale io ~ I/Kq. Inserting the expression for the energy flux @ inside relation 
(O), one obtains the 2D EMHD counterpart of Kolmogorov's four-fifths law il|J 



Q3(^) ^ ei (12) 

where 

Q^ii) = -\dl S,{i) - \dl V,{i) + 

I 1 2 1 2 (13) 

Relation (0) relies only on the aforementioned hypothesis here recalled: homogeneity, 
isotropy, and the existence of an inertial range of wavenumbers in which the energy flux 
is constant, with a value tending to a finite positive limit for vanishing viscosity. The most 
remarkable aspect of relation ([T^JT3| ) is that it does not only provide a linear scaling for 
the third order structure function Q3{i) within the inertial range of lengthscales, but it also 
prescribes the value of the numerical coefficient appearing in front of the scaling relation. 
Moreover, it is valid for any value of de, meanwhile no power law scaling relation is expected 
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to hold for, say, the second order structure functions, apart from hmiting cases such as 
4 ~ 1 and 4 < 1 dig. 

To check the vahdity of EMHD Kolmogorov law (p!^), the equations (|I]) are solved in a 
square box of size 27r x 27i imposing periodic boundary conditions, by means of a standard 
pseudospectral method with resolution N x N. Hyper dissipation with g = 4 is employed in 
order to achieve a larger extent of the inertial range, which in physical space is known to 
be much narrower than in spectral space 0. Hyperviscosity is set to z/4 = 10~^^, 10~^^ for 
= 512 simulation and z/4 = 10"^^ for = 1024. 

The initial conditions are v{k) ~ fc^ exp(— A;^/2/cq) with random phases, ko = 1 and total 
energy of order unity in both resolution simulations. After a transient of a few large eddy 
turnover times, when the energy initially contained at the lowest wavenumbers starts to 
cascade down to small scales, the energy dissipation reaches a maximum value and then a 
self similar stage of decay sets in [0. The energy flux is approximately constant throughout 
the inertial range of wavenumbers (see Fig.l), and, in agreement with the assumption (^, 
its value appears to be asymptotically independent of viscosity. The structure function 
Qsli) is computed during the self similar stage of decay. The results are obtained after 
averaging over a short time in order to get better statistics at small scales. As shown in 
Fig. 2 the compensated structure function Qs{i)/{ei) approaches unity, as prescribed by the 
relation (0) in an interval delimited from below by the "Taylor scale" A and above by the 
energy containing scale Iq. By lowering the viscosity (crosses, A^ = 1024) the scaling range 
extends to smaller scales over almost one decade. As previously remarked, the width of the 
inertial range is actually diminished by the fact that, for a decaying flow, the energy flux 
in wavenumber space is not constant except asymptotically (see Fig.l). This results are an 
evident numerical confirmation of the validity of the Kolmogorov-type relation (|12|). This 
kind of assessment is important since it lies at the foundations of the statistical study of 
turbulence. 

Introducing the further hypothesis of statistical self-similarity one can infer the following 
scaling behavior for the velocity and for the magnetic field differences in the asymptotic case 
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^ < 4 , (14) 
which leads to the small scales energy spectrum, dominated by kinetic energy, 

E{k) = CKS^/^k-''/^ kde^l. (15) 

On the other hand, for scales larger than c/g, the expected self-similar scaling is 

^ » 4 (16) 

5B{1) OC £i/3£2/3 
leading to the energy spectrum, dominated by magnetic energy, 

E{k) = C;,£2/3^-7/3 A;4 < 1 . (17) 

The slopes of the computed spectra, as shown in Fig. 3, are close to the estimates (|15[) and 
(P^) over a wide range of wavenumbers. 

The main results of this work are the derivation of a Kolmogorov-type relation for 2D 
EMHD decaying turbulence, and its numerical confirmation. Since the Kolmogorov's law 
stands as a corner stone in the study of the statistical features of turbulence, these results 
form the basis of further analysis, starting from the issue of intermittency. 

The authors wish to acknowledge support and hospitality by the Istituto di Cosmoge- 
ofisica CNR, Torino. The calculations were partially performed with computer facilities of 
INFN, Sezione di Torino. 
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FIGURES 

FIG. 1. Energy flux in wavenumber space normalized to energy dissipation, (a) = 0.3. 
Continuous line: = 512; 1/4 = IQ-^^ Dotted line: = 1024; 1/4 = 10~^^. (b) 4 = 0.02. 
Continuous line: N = 512; z^4 = 10~^^. 

FIG. 2. Compensated structure function Q3{i)/{£i)- (a) = 0.3; Diamonds: N = 512; 
1/4 = 10-12. Crosses: N = 1024; 1/4 = lO'^^ ^ q Diamonds: N = 512; U4 = lO^^^. 

FIG. 3. Energy spectrum, (a) 4 = 0.3; Diamonds: N = 512; 1/4 = lO'^^^ Crosses: N = 1024; 
1/4 = 10^1^. The continuous line is the Kolmogorov spectrum ( [T5| ) with Ck = 2.0. (b) = 0.02; 
Diamonds: N = 512; 1/4 = lO^^"^. The continuous line is the spectrum (^) with C'j^- = 8.0 
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